Small-x physics and heavy quark 
photoproduction in the semihard approach at 

HERA. 



V.A.Saleev 

Samara State University, Samara 443011, Russia 

N.P.Zotov 

Nuclear Physics Institute, Moscow State University, Moscow 119899, Russia 



Abstract 

Processes of hevy quark photoproduction at HERA energies and be- 
yond are investigated using the semihard {k± factorization) approach. 
The virtuality and longitudinal polarization of gluons in the photon - 
gluon subprocess as well as the saturation effects in the gluon distribu- 
tion function at small x have been taken into account. The total cross 
sections, rapidity and p± distributions of the charm and beauty quark 
photoproduction have been calculated. We obtained the some differences 
between the predictions of the standard parton model and the semihard 
approach used here. 



1 Introduction 



The heavy quark production at HERA is very intresting and important subject of study 
^, H Because it is dominated by photon-gluon fusion subprocess (Fig.l) one can study 
the gluon distribution functions G{x, Q^) in small x region (roughly at x > 10"''). Secondly 
this last issue is important for physics at future colliders (such as LHC): many processes at 
these colliders will be determined by small-x gluon distributions. 




Fig. 1 

The heavy quark production processes are so-called semihard processes ||^ at HERA en- 
ergies and beyond. By definition in these processes a hard scattering scale Q (or heavy 
quark mass M) is large as compared to the Aqcd parameter, but Q is much less than the 
total center of mass energies: Aqcd ^ Q <C y/s. Last condition implies that the processes 
occur in small x region: x ~ M'^/s <C 1. In such case parturbative QCD expansion have 
large cofScients 0(ln"(s/M^)) ~ 0(ln"(l/a;)) besides the usual renormalization group ones 
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which are 0{\n^ (Q"^ / h?)) . So in perturbative QCD the heavy quark photoproduction cross 
section as result of photon-gluon fusion processe has the form a^g = a°g + ctgal^g + 
where al^g was calculated by Ellis and Nason j^]. The photon - gluon fusion cross sec- 
tion in low order decreases vs s at s — > oo : a°g ^ M^/s ln(s/Q^). But in the same limit 
a}^g — > Const, because of that the heavy quark photoproduction cross section is dominated 
by the contribution of the gluon exchange in the t - channel. That results breakdown the 
standard perturbative QCD expansion and the problem of summing up all contributions of 
the order (a^ ln(QVA2))", (a^ ln(l/x))" and (a^ ln(QVA2) ln(l/a;)" in perturbative QCD 
appears in calculation of a^g. 

It is known that summing up the terms of the order [as ln((3^/M^))" in leading logarithm 
approximation (LLA) of perturbative QCD leads to the linear DGLAP evolution equation 
for deep inelastic structure function . Resummation of the large contributions of the order 
of {us ln(l/a:))" leads to the BFKL evolution equation |^ and its solution gives the drastical 
increase of the gluon distribution: xG{x,Q^) ^ x^'^'^ujo = (4Afc ln2/7r)as(Q§). The speed 
of the growth of the parton density as a; ^ makes the parton-parton interactions very 
important, which in turn makes the QCD so-called GLR evolution equation essentially 
nonlinear . 

The growth of the parton density at x — > and interaction between partons induces sub- 
stantial screening (shadowing) corrections which restore the unitarity constrains for deep 
inelastic structure functions (in particular for a gluon distribution) in small x region | |lO| . 
These facts break the assumption of the standard parton model (SPM) about x and trans- 
verse momentum factorization for a parton distribution functions. We should deal with the 
transverse momentum factorization (k± factorization) theory |ll|, Making allowance 
for screening corrections at small x we have so-called semihard approach ^ , which take 
into account the virtuality of gluons, their transverse momenta as well as the longitudional 
polarizations of initial gluons. 

In semihard approach j^, |lOj screening corrections stop the growth of the gluon distribution 
at a; — !■ 0. This effect was interpreted as saturation of the parton (gluon) density: gluon 
distribution function xG{x,Q'^) becames proportional to Q^R^ at < q^x) and a ~ 
-^xG{x, Q^) ~ R^. The parameters R and qgix) can be considered as new phenomenological 
parameters: the R is related to size of hadron or of black spots in a hadron and the parameter 
Qq is a typical transverse momentum of partons in the parton cascade of a hadron, which 
leads to natural infrared cut-off in semihard processes. The parameter Qq increases with 
growth of s H^, |l3|. So, because of the high parton density in small x region a standard 
methods of perturbtive QCD can not be used 10, ^ In this point the semihard 
approach ||^, |l^ differs from the works |l5| where after choice of structure function 

at starting value Qq the QCD evolution are calculated using certain equations. 

In the region, where the transverse mass of heavy quark M± << y^, one need take into 
account the dependence of the photon-gluon fusion cross section on the virtuality and polar- 
izations of initial gluon. Thus in semihard approach the matrix elements of this subprocess 
differs from ones of SPM (see for example |l2|). 

We would like to remark that as corrections in the semihard approach look quite different 
in comparison with usual ^ and k± factorization n, |l^ ones. This corrections have 
been taken into account using new gluon density ^{x, gj) (see below) which was calculated 
by Gribov, Levin and Ryskin |^ in LLA where all large logs (of the types ln{l/x) and 
ln{Q'^ / X^)) were summed up. The another corrections which give contributions to if —factor 
have been taken into account in the normalization of function xG{x,Q^) (see below also) 

0- 
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As far as the SPM calculations of the next-to-leading (NLO) cross section to the photo- 
production of heavy quarks that the review of its may found in the papers The 
results 0] have been confirmed by Smith and Van Neerven The futher results for 

the electro- and photoproduction of heavy quarks was obtained in Refs. |T^. Authors 
in ]l6| notes that for beauty quark production the NLO corrections are large and various 
estimates of corrections lead to theoretical uncertainties of the order of a factor of 2 to 3. 
As was estimated in the total cross section for beauty quark production at HERA will 
be only a few tens of percent large than the one- loop results Q . 

For charm quark a theoretical uncertainties are even the higher. It is known that there is 
also the strong mass dependence of the results of calculations for charm quark production. 
Since the mass of the charm quark is small for perturbative QCD calculations the resumme 
procedure ||ll|] is not applicable for charm quark production . 

In Refs. ^ we used the semihard approach in order to calculate the total and differ- 
ential cross sections of the heavy quarkonium, J/vl/ and T, photoproduction. We obtained 
the remarkable difference between the predictions of the semihard approach and the SPM 
especially for p±- and ^-distributions of the J/'if mesons at HERA. 

In present paper we investigate the open heavy quark production processes in the semihard 
approach, which was used early in Rcf. [ic| ] for calculation of heavy quark production rates 
at hadron colliders and for prediction of J/^* andT photoproduction cross sections at high 
energies in Refs. p^ . 



2 QCD Cross Section for Heavy Quark Electroproduction 

We calculate the total and differential cross sections (the p± and rapidity distributions) 
of charm and beauty quark photoproduction via the photon-gluon fusion QCD subprocess 
(Fig.l) in the framework of the semihard (fcj_ factorization) approach First of all 

we take into account the transverse momentum of gluon q2±-, its virtuality g| — — (zfj^ 
and the alignment of the polarization vectors along the transverse momentum such as 
- 92±m/ I 9*2^ I ©111. 

e{Pi) 



P{P2) 

Fig. 2 

Let us define Sudakov variables of the process ep QQX (Fig. 2): 

Pi = ctiPi + /3iP2 + pi±, P2 ^2?! + P2P2 + P2^, 

qi = XiPi+qi±, q2=X2P2 + q2±, (1) 
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where 



2 2 71 /r2 2 2 2 2 



Pi and p2 are 4-monienta of the heavy quarks, qi is 4-momentum of the photon, q2 is 4- 
niomentum of the gluon, pij., P2±, qi±^ 92_l are transverse 4- momenta of these ones. In the 
center of mass frame of cohiding particles we can write Pi = {E, 0, 0, E), P2 = {E, 0, 0, —E), 
where E ~ \/s/2, Pi = P| = and (P1P2) = s/2. Sudakov variables are expressed as 
follows: 

a\ = ^^exp(2/J, a2 = ^^exp(y2), 

01 = ^^exp{~yl), (32 ^ ^^^exp{-y;,), (2) 

where Mf 2± = ^'^^ + Pi 2±j 2/* 2 are rapidities of heavy quarks, M is heavy quark mass. 
From conservation laws we can easy obtain following conditions: 

qi± + q2± = Pi± + P2±, xi=ai+a2, X2^l3i+P2- (3) 
The differential cross section of heavy quarks electroproduction has form: 

(ep ^ QQX) = / dyldy; ^ ^ |MP<i>.(x, g^)<l> (..,,2^) 



(Ppi± J IT TT 167r^(a;ia;2s)^ 

For photoproduction process it reads: 



(Ppi±_ J TT 16TT^{sX2)^a2 

We use generalized gluon structure function of a proton ^p{x2,q2±) which is obtained in 
semihard approach. When integrated over transverse momentum q2i_ (92^ = (0,(f2±,0)) 
of gluon up to some limit it becomes the usual structure function giving the gluon 
momentum fraction distribution at scale Q^: 

J %ix, ql^)df2i_ = xGpix, Q^). (6) 
We use in our calculation following phenomenological parameterization px| ]: 

%{^, <ll) = - xfUix, gi), (7) 

where 

fp = 1, ?! < <ll{x) 

U - (^)^ qi>ql{:<^\ (8) 

and ql{x) = + exp(3.56 v'ln(xo/a;)), Ql = 2GeV^ A = 56 MeV, xq = 1/3. The 
normalization factor C ~ 0.97 mb of the structure function <i>p(a;, (?f^) was obtained in [|lOj 
where 66-pair production at Tevatron energy was described. 

The $e(a^i; ^ij,) is well known virtual photon spectrum in Weizsacker- Williams approxima- 
tion psi before the integration over (7^: 
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The effective gluon distribution xG{x, Q^) obtained from (7) increases at not very small x 
(0.01 < a:; < 0.15) as x~'^° , where ojq — 0.5 corresponds to the BFKL Pomeron singularity 
[^.This increases continuously wp to x = xq, where xq is a solution of the equation (/qC^^o) — 
Q^. In the region x < xq there is the saturation of the gluon distribution: xG{x, Q^) ~ CQ^. 

The square of matrix element of partonic subprocess 7*5* — > QQ can be written as follows: 



\M\ — Idn eQasa{xiX2s) 



For real photon and off-shell gluon it reads: 



' + (10) 



(i-M2)(u-M2) 

where a2 = 1 — ai and s, i, u are usual Mandelstam variables of partonic subprocess 

S = bl+P2)' = (gi+92)', t=(pi-9i)2 = (p2-92)', (12) 

u = (Pl-g2)^ = (p2-9l)^ s + < + 'u = 2M2 + g2_^+g2_^. 



3 Discussion of the Results 



The results of our calculations for the total cross sections of c- and 6-quark photoproduction 
are shown in Fig. 3 (solid curves). Dashed curves correspond to the SPM predictions with 
the GRV parametrization of the gluon distribution |Q . The results of calculations for the 
SPM are shown without K-factor, wich have the typical value K = 2 for hard scattering 
processes. (In the semihard approach K-factor is absent [|l0[ ^). The K-factor in the 
SPM may change the relation between the results of calculations in the semihard approach 
and SPM. But independently from it we see that the saturation effects more clearly are 
pronounced for charm quark photoproduction (at y/Sjp ~ 500 GeV). In any case (with 
or without K-factor in SPM) the cross section for beauty quark photoproduction in the 
semihard approach is the higher than one in the SPM at y/s^ > 200 GeV. 

The p± distributions for c- and 6-quark photoproduction in the semihard approach (solid 
curves) and in the SPM (dashed curves) are shown in Fig. 4. The curves are obtained in 
the semihard approach for charm quark photoproduction show the saturation effects in low 
P± region {pj_ < 2 Gev/c). In middle p± region (2 GeV/c < p± < 15 GeV/c) the heavy 
quark photoproduction pj_ distributions in the semihard approach are the higher ones of 
the SPM (with GRV paramatrization of gluon distribution). At high > 15 GeV/c we 
have contrary relation between p±^ distributions in the semihard approach and SPM. 



Thus the SPM leads to over-estimatedcross section in the low p± region and under- estimated 
one in the middle p± region as was noted as early as in Refs. |10, 15| (see also ||2^). This 
behavior of p± distributions in the k± factorization approach is result from the off mass shell 
subprocess cross section as well as the saturation effects of gluon structure function in 
semihard approach [hol. 
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Fig. 5 show the comparison of the; rc;siihs for heavy quark rapidity distribution (in the 
photon-proton center of mass frame) in the different models: sohd curve shows the y dis- 
tribution in the semihard approach, dashed curve shows one in SPM. The discussed above 
effects arc sufficiently large near the kinematic boundaries, i.e. at big value of We see 
that the difference between solid and dushed curvers can't be degrade at all y* via change 
of normalization of SPM prediction. 



4 Conclusions 



We shown that the semihard approach leads to the saturation effects for the total cross 
section of charm quark photoproduction at available energies and predicts the remarkable 
difference for rapidity and transverse momentum distributions of charm and beauty quark 
photoproduction, which can be study already at HERA ep collider. 
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Figure captions 

1. QCD diagrams for open heavy quark photoproduction subprocesses 

2. Diagram for heavy quark elecrtoproduction 

3. The total cross section for open charm and beauty quark photoproduction: solid curve 
- the semihard approach, dashed curve - the SPM 

4. The p± distribution for charm and beauty quark photoproduction at = 200 



GeV: curves as in Fig. 3 
5. The y* distribution for charm and beauty quark photoproduction: at y^s-yp = 200 



GeV: curves as in Fig. 3 
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